Unexpected Nonlinear Dynamics in NbN Superconducting Microwave Resonators 
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In this work we characterize the unusual nonhnear dynamics of the resonance response, exhib- 
ited by our NbN superconducting microwave resonators, using different operating conditions. The 
nonhnear dynamics, occurring at relatively low input powers (2-4 orders of magnitude lower than 
Nb), and which include among others, bifurcations in the resonance curve, hysteresis loops and reso- 
nance frequency shift, are measured herein using varying temperature, applied magnetic field, white 
noise and rapid frequency sweeps. Based on these measurement results, we consider a hypothesis 
according to which Josephson junctions forming weak links at the boundaries of the NbN grains 
are responsible for the observed behavior, and we show that most of the experimental results are 
qualitatively consistent with such hypothesis. 

PACS numbers: 85.25.-j, 74.50.+r, 47.20.Ky, 05.45.Xt 



I. INTRODUCTION 

Understanding the underlying mechanisms that cause 
and manifest nonhnear effects in superconductors has a 
significant importance both theoreticaUy and practically. 
On the one hand, knowing how to minimize and control 
nonlinear effects, being the major cause for the degra- 
dation of performance of RF devices employing super- 
conductors 0, would have a great technological impact. 
It will enable the design and fabrication of high perfor- 
mance superconducting RF devices with extended linear 
behavior, and higher power densities than afforded nowa- 
days 0,13; On the other hand, identifying the origins of 
the nonlinearities in superconductors, would lead to a 
better understanding of the superconducting phenomena 
and its behavior in the microwave regime. 

In spite of the intensive study of nonlinearities in su- 
perconductors in the past decades, we still lack a coherent 
picture regarding these effects and their origins Q . This 
is partly because the nonlinear mechanisms in supercon- 
ductors are various and many times act concurrently on 
the superconducting sample, hence making the identifi- 
cation process of the dominant factor mainly indirect and 
based on eliminations . 

Nonlinear mechanisms in superconductors, which are 
usually divided into intrinsic and extrinsic origins, in- 
clude among others, Meissner effect Ig, pair-breaking, 
global and local heating effects 0, 137 rf and dc vor- 
tex penetration and motion |9|| , d efect points, damaged 
edges ^3 1 substrate material and weak links [l^ . 
Whereas weak links is a collective term representing var- 
ious material defects located inside the superconductor 
such as weak superconducting points switching to nor- 
mal state under low current density, Josephson junc- 
tions forming inside the superconductor structure, grain- 
boundaries, voids, insulating oxides, insulating planes. 
These defects and impurities generally affect the con- 
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duction properties of the superconductor and as a result 
cause extrinsic nonlinear effects. 

Nonlinear effects in superconductors in general and 
in NbN in particular have been reported and analyzed 
by several research groups. Duffing like nonlinear! ty for 
example was observed in different superconducting res- 
onators employing different geometries and materials. It 
was observed in a HTS parallel plate resonator j\, in a 
Nb microstrip resonator 13],inaNb and NbN stripline 
resonators [3] , in a YBCO coplanar- waveguide resonator 
[l^ . in a YBCO thin film dielectric cavity and also 
in a suspended HTS thin film resonator [l5j. Other 
nonlinearities including notches, anomalies developing at 
the resonance lineshape and frequency hysteresis can be 
found in 

In this paper we report the observation of unex- 
pected nonlinear dynamics measured in NbN supercon- 
ducting microwave resonators. Among the observed ef- 
fects, asymmetric resonances, multiple bifurcations in the 
resonance lineshape, hysteretic behavior in the vicinity of 
the bifurcations, hysteresis loops changing direction, res- 
onance frequency shift as the input power is increased 
and nonlinear coupling. Some of these nonlinear ef- 
fects were introduced in a previous publication [T^ , thus 
this paper will concentrate on presenting a different set 
of measurements applied to these nonlinear resonators, 
which provides a better understanding of the underlying 
physical mechanism causing these effects. For this pur- 
pose, we have measured the nonlinear superconducting 
resonators under different operating conditions, such as 
added white noise (to the main signal), fast frequency 
sweep using frequency modulation (FM), applied mag- 
netic field, and varying temperature. Under each one 
of these operating conditions we observe interesting non- 
linear dynamics, which qualitatively agree with our hy- 
pothesis of microscopic Josephson junctions forming at 
the boundaries of the NbN columnar structure. 

The remainder of this paper is organized as follows, 
the fabrication process of the NbN superconducting res- 
onators is described briefly in Sec. II. A short summary 
of previous results is brought in Sec. III. The nonlin- 
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FIG. 1: (a) Schematic cross section of the stripline resonator 
used, consisting of five layers: two superconducting ground 
planes, two sapphire substrates and a NbN film in the middle 
deposited on one of the sapphires, (b) Top view of the three 
resonator layouts (Bl, B2, B3) which were deposited as the 
middle layer. 



ear dynamics of these microwave resonators measured 
under various operating conditions are reviewed in Sec. 

IV. Whereas possible underlying physical mechanisms re- 
sponsible for the measured effects are discussed in Sec. 

V, which also concludes this paper. 
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FIG. 2: Sii amplitude measurement of Bl resonator at its 
first mode ~ 2.59GHz. The measured resonance lineshapes 
are asymmetrical and contain two abrupt bifurcations at the 
sides of the resonance. Moreover the bifurcation frequencies 
shift outwards as the input power is increased. The measured 
resonance lineshapes were shifted vertically by a constant off- 
set for clarity. The layout of Bl resonator is depicted in the 
inset. 



III. SUMMARY OF PREVIOUS RESULTS 



II. FABRICATION PROCESS 



The measurement results presented in this paper 
belong to three nonlinear NbN superconducting mi- 
crowave resonators. The resonators were fabricated us- 
ing stripline geometry, consisting of two superconducting 
groimd planes, two sapphire substrates, and a deposited 
strip in the middle (the deposition was done on one of 
the sapphire substrates). Fig. shows a schematic di- 
agram illustrating stripline geometry and a top view of 
the three resonators layouts. For convenience we will 
be referring to the three resonators in the text by the 
shortened names Bl, B2 and B3 as defined in Fig. ^ 
The sapphire substrates dimensions used were 34mm X 
30mm X 1mm, whereas the coupling gap between the res- 
onators and their feedline was set to 0.4mm in Bl and B3 
cases and 0.5mm in B2 resonator. The resonators were 
dc-magnetron sputtered in a mixed Ar/N2 atmosphere, 
near room temperature. The patterning was done us- 
ing standard UV photolithography process, whereas the 
NbN etching was performed by Ar ion-milling. The sput- 
tering parameters, design consideration as well as other 
fabrication details can be found elsewhere |0 . The crit- 
ical temperature Tc of Bl, B2 and B3 resonators were 
relatively low and equal to 10. 7K, 6.8K and 8.9K respec- 
tively. The thickness of the NbN resonators were 2200 A 
in Bl, 3000 A in B2, and 2000 A in B3 resonator. 



As we have shown comprehensively in a previous pub- 
lication two of the main nonlinear features observed 
in our NbN nonlinear superconducting resonators are bi- 
furcations that appear in the resonance curves at rela- 
tively low input powers and hysteresis loops forming in 
the vicinity of the bifurcations. In Fig. |2] we present 
a 5*11 parameter measurement of Bl first mode reso- 
nance featuring bifurcations at both sides of the reso- 
nance curve. 

Furthermore, the frequencies at which the bifurcations 
occur shift outwards from the center frequency as the in- 
put power is increased. The measurement was performed, 
at liquid helium temperature 4.2K, using vector network 
analyzer. The measured curves in the figure were shifted 
vertically by a constant offset for clarity. In contrast to 
the nonlinear dynamics presented in Fig. 12 we show in 
Fig. 131 a resonance response measured at 4.2K, exhibit- 
ing Duffing nonlinearity of the kind generally reported 
in the literature [3, ll4L Il5| . This nonlinearity which can 
be explained in terms of resistance change AR and ki- 
netic inductance change ALk HOl was measured at 
the first resonance frequency of a 2200 A thickness Nb 
resonator employing B2 layout geometry (Tc — 8.9K). 
The differences between the two nonlinear dynamics are 
obvious, but nevertheless note the difference in the or- 
der of magnitude of the input powers at which these two 
nonlinearities occur, ^ —28 dBm in the NbN case vs. 
~ 10 dBm in the Nb resonator. In Fig. 01 we present 
a Sii measurement of B2 resonator at its second mode, 
measured while sweeping the frequency in both direc- 
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FIG. 3: DufRng nonlinearity exhibited by a Nb stripline res- 
onator employing B2 layout at its first mode. The differ- 
ent 511 amplitude plots correspond to different input powers, 
ranging from —15 dBm to 15 dBm in steps of 2 dBm. As the 
input power is increased the resonance becomes asymmetri- 
cal and infinite slope develops at the left side of the resonance 
curve. The plots were offset in the vertical direction for clar- 
ity. 



tions. The red line represents a forward frequency sweep 
whereas the blue Une represents a backward frequency 
sweep. The graphs corresponding for different input pow- 
ers were also offset in the vertical direction for clarity. At 
—8.04 dBm the resonance is linear and shows no hystere- 
sis behavior. As the input power is increased by a power 
step of 0.01 dBm, to —8.03 dBm the resonance response 
lineshape is changed dramatically, featuring two bifurca- 
tions at both sides of the resonance curve and hysteresis 
loops forming in the vicinity of the bifurcations (as can 
be clearly seen in the figure) . Moreover this measurement 
exhibits another interesting nonlinear feature. The right 
side hysteresis loop changes direction as the input power 
is increased. Between —8.03 dBm and —7.98 dBm the 
hysteresis loop circulates counter clockwise. At —7.98 
dBm the hysteresis loop vanishes. Whereas as the input 
power is increased further, the loop circulates clockwise. 
Furthermore at —6.35 dBm (not shown here) the hys- 
teresis loop vanishes again and changes its direction one 
more time, (it starts circulating counter clockwise as the 
power is increased). 

In Fig. |S1 we show yet another nonlinear feature ex- 
hibited by these nonlinear resonators, namely multiple 
bifurcations in the resonance curve. The figure plots the 
resonance response of B3 resonator at its first resonance 
frequency, corresponding to 1.49 dBm input power. The 
measurement which was obtained by sweeping the fre- 
quency axis in the forward and backward directions, fea- 
tures three bifurcations and hysteresis loops within the 
resonance lineshape in each direction. 

Other experimental data featuring these effects and 
similar nonlinearities observed in Bl, B2 and B3 res- 
onators can be found in Ref. 1 181 . Among the nonlinear- 
ities not brought here or in |l8l |. one can name, nonlinear 
coupling '2ll|, and intermodulation gain [2^ . 
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FIG. 4: Frequency sweep measurement of B2 resonator at 
~ 4.395GHz performed in both frequency directions. The 
plots exhibit hysteresis loops forming at the vicinity of the 
bifurcations, as well as hysteresis loops changing direction as 
the input power is increased. The red line represents a forward 
sweep, whereas the blue line represents a backward sweep. 
The resonance lineshapes were shifted vertically by a constant 
offset for clarity. 



IV. MEASUREMENT RESULTS 

In the following subsections we present measurement 
results of Bl and B2 resonators, under different operat- 
ing conditions. In subsection A, the impact of white noise 
power on the bifurcation feature is briefly examined. The 
dependence of the nonlinear dynamics on temperature 
varying and applied magnetic field are investigated in 
subsections B and C respectively. In subsection D, heat- 
ing effects are considered through frequency sweep time 
analysis and comparisons to heating time scales previ- 
ously reported in the literature. Evidence of the colum- 
nar structure of our NbN films is presented in subsection 
E. Possible interpretations of these observations are dis- 
cussed in Sec. V. 



A. WHITE NOISE EFFECT 

Bifurcations in the resonance response of nonlinear os- 
cillators is usually described in terms of metastable states 
and dynamic transition between basins of attraction of 
the oscillator [23] . Thus in order to examine qualitatively 
the metastability of the resonance bifurcations of these 
nonlinear resonators, we have applied a constant white 
noise power to the resonator, several orders of magni- 
tude lower than the main signal power, using the setup 
shown in Fig. |B1 The white noise level applied was —58 
dBm/Hz. It was generated by amplifying the thermal 
noise of a room temperature 50 f2 load using an am- 
plifying stage. The generated noise was added to the 
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FIG. 5: Sii parameter amplitude measurement of the first 
resonance of B3 resonator, (shown in the inset), measured at 
input power of 1.49 dBm. The measurement was done using 
a network analyzer employing 4000 measurement points in 
each direction. The red line represents a forward frequency 
scan whereas the blue line represents a backward scan. The 
plot shows clearly three bifurcations within the resonance line- 
shape in each direction, as indicated by small circles. 



transmitted power of the network analyzer via a power 
combiner. The power reflections were redirected by a 
circulator on their way back, and were measured at the 
second port of the network analyzer, thus measuring 5*21 
parameter. The effect of the —58 dBm/Hz white noise 
power on Bl first mode bifurcations, is shown in Fig. 
{7\ (a), whereas for comparison in Fig. [T] (b) we show 
nearly noiseless —80 dBm/Hz) resonance curve mea- 
surements performed using the simple setup depicted in 
Fig. El without the amplifier and the combiner stage. 
The two measurements were carried out within the same 
input power range (—23.9 dBm through —20 dBm). 

By comparing between the two measurement results, 
one can observe the following, the two fold bifurcations 
in Fig. [3 (b) form a hysteresis loop at each side of the 
resonance curve. In Fig. |7| (a) in contrast, as a result of 
the added noise, the bifurcations at the left side, repre- 
sented by the thick colored lines on the graphs, become 
frequent and bidirectional, whereas the hysteresis loops, 
at the right side, vanish. 

The transition rate T (/) between the oscillator basins 
of attraction (as a function the oscillator frequency /), 
can be generally estimated by the expression 23] T (/) — 
roexp{-EA{f) /kBTeff) , where E'a (/) is the quasi- 
activation energy of the oscillator, Te// is proportional 
to the noise power, is Boltzmann's constant, / is the 
oscillator frequency, whereas Fq is related to Kramers 
low-dissipation form [2^ and it is given approximately by 
fo/Q, where fo is the natural resonance frequency, and 
Q is the quality factor of the oscillator. Based on the 
results presented in Fig. [7| one can evaluate the foUow- 
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FIG. 6: Schematic diagram of the experimental setup used 
to measure the nonlinear resonance 2.585 GHz of Bl using 
frequency sweep mode while applying —58 dBm white noise. 
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FIG. 7: Frequency sweep measurement of Bl resonator at 
2.585 GHz performed in both directions (a) while applying 
white noise of —58 dBm/Hz (b) without applying external 
noise (~ —80 dBm/Hz). The red line represents a forward 
sweep, whereas the blue line represents a backward sweep. 
The measured resonance curves were shifted vertically by a 
constant offset for clarity. 



ing parameters fo ~ 2.585GHz, Q - 1620, Fq - 1.6MHz, 
Te// ~ lO^'^K. Whereas Ea the quasi-activation energy 
of the oscillator may be roughly estimated for the left and 
right sides of the resonance response separately. For the 
right side of the resonance curve, the quenching of the 
hysteresis as a result of the white noise, suggests that 
the applied noise power —58 dBm/Hz is within the high 
noise limit, where the noise power is higher than the en- 
ergy barrier associated with this bifurcation. Whereas 
for the left side, the hysteresis loop and the noise-induced 
transitions, indicate that the noise power —58 dBm/Hz 
is still in the moderate noise limit, where the noise power 
is less than the quasi-activation energy associated with 
this bifurcation. 



B. TEMPERATURE DEPENDENCE 

The physical properties of superconductors such as 
London penetration depth, surface resistance and kinetic 
inductance are strongly dependent on temperature |2^. 
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Whereas temperature varying measurements are gener- 
ally used in the characterization process of supercon- 
ductors and in the determination of London penetration 
depth [2^1 , we vary the temperature here in order to ex- 
amine its effect on the nonlinear dynamics of the res- 
onators. The measurements were performed at temper- 
ature ranges well below Tc of the resonators, thus the 
exhibited nonlinear dynamics presented here are features 
of the superconducting phase only. In Fig. |S1 subplot (a) 
we show a iSii parameter amplitude measurement ver- 
sus frequency of Bl resonator attained at its third res- 
onance mode 0. The input power was set to —16.9 
dBm, a power level at which the bifurcation is already 
present in the resonance curve, afterwards the temper- 
ature was increased from 5.4K through 8K in steps of 
O.IK. The nonlinear dynamics of the resonance curve as 
the temperature is increased exhibit three main nonlin- 
ear features, the resonance frequency shifts towards lower 
frequencies, the bifurcation step decreases and the bifur- 
cation frequency becomes lower. Whereas the first non- 
linear feature can be explained in terms of increase in 
the resonator inductance per unit length L due to tem- 
perature rise [2^, the second and third features are less 
straightforward, and are more likely related to the effect 
of temperature on the underlying nonlinear mechanism 
responsible for the observed effects. Likewise in subplot 
(b), B2 resonator at its second mode exhibits a similar 
behavior, though the input power applied in this case is 
— 10 dBm, which is lower than the bifurcation threshold 
of this resonance ( ~ — 9 dBm). Similarly the resonance 
frequency shifts gradually towards lower frequencies as 
the ambient temperature is varied from 5.4K through 
5.9K in steps of O.OIK. Moreover at 5.8K, the resonator 
even achieves critical coupling condition (^n ~ 0) via 
temperature varying only. 
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FIG. 8: Subplot (a) exhibits the nonlinear resonance fre- 
quency response of Bl measured under constant input power 
of —16.9 dBm while increasing the temperature from 5.4K 
through 8K in steps of O.IK. In addition to the gradual reso- 
nance frequency shift, one can notice the resonance evolution 
as the temperature increases. Subplot (b) shows B2 second 
mode measured at a constant input power of —10 dBm while 
increasing the temperature from 5.4K through 5.9K in steps of 
O.OIK. In addition to the gradual resonance frequency shift, 
one can notice a temperature induced critical coupling occur- 
ring at 5.8K. The resonance curves were shifted vertically by 
a constant offset for clarity. 



the resonance response while increasing the input power. 
The dynamic behavior measured using this method, as 
the input power is increased, undergoes through the fol- 
lowing sequential phases, symmetrical and Lorentzian 
curves, resonance curves containing an upward bifurca- 
tion, a curve without bifurcation, and finally resonance 
curves containing a downward bifurcation. 



In order to better examine the extrinsic origin of the 
observed nonlinearities, we have investigated the depen- 
dence of the nonlinear resonances on magnetic field. The 
magnetic field dependence was measured by applying two 
methods, one by setting a constant input power and vary- 
ing the magnetic field, and two by setting a constant mag- 
netic field and varying the input power. The magnetic 
field applied was perpendicular to the resonator plane. 
The magnetic field measurements exhibiting B2 second 
resonance mode are summarized in Figs. |31 IIOI In or- 
der to obtain the results of Fig. |^ we have applied a 
constant input power of —5 dBm to the resonator, and 
measured its Sn response while increasing the applied 
magnetic field by small steps. Above some low magnetic 
field threshold of II.SmT, the left side bifurcation van- 
ishes, thus suggesting that the physical mechanism caus- 
ing the bifurcations is sensitive to low magnetic field. In 
Fig. ^|we applied the second method, where we set the 
magnetic field to a constant level of 0.09T and measured 



D. FREQUENCY SWEEP TIME ANALYSIS 

Resistive losses and heating effects are typically char- 
acterized by relatively long time scales |3;ll3- In attempt 
to consider whether such effects are responsible for the 
observed nonlinearities in general and for the bifurcations 
in particular, we have run frequency sweep time analysis 
using the experimental setup depicted in Fig. We 
have controlled the frequency sweep cycle of an Anritsu 
model signal generator via FM modulation. The FM 
modulation was obtained by feeding the Anritsu with a 
saw tooth waveform generated by another signal gener- 
ator having 1// sweep time cycle. The reflected power 
from the resonator was redirected using a circulator and 
measured by a power diode and oscilloscope. The left 
and right hand bifurcations of B2 ~ 4.39GIIz resonance 
were measured using this setup, while applying increasing 
FM modulation frequencies: 1, 10, 20, 30, 40, 50, 100, 150 
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FIG. 9: Increasing the magnetic field gradually from zero 
causes the bifurcation in the B2 resonance band to disappear 
at relatively low value of ll.SmT, while applying a constant 
input power level of —5 dBm. This bifurcation vanishing indi- 
cates that the bifurcation mechanism is sensitive to magnetic 
field. The measured curves were shifted vertically by a con- 
stant offset for clarity. 
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FIG. 10: B2 nonlinear resonance response measured under 
a constant magnetic field of 0.09T, while increasing the in- 
put power. The figure exhibits additional dependencies of 
the resonance response curves on applied magnetic field. The 
resonance which starts as a Lorentzian develops into a res- 
onance curve having an upward bifurcation as the power is 
increased, afterwards into a curve with no bifurcation, fol- 
lowed by a resonance curve having a downward bifurcation as 
the input power is increased further. The measured curves 
were shifted vertically by a constant offset for clarity. 



and 200kHz. In Fig. ^] we present a measurement 
result obtained at 50kHz FM modulation, or alterna- 
tively Tsweep of 20fis. The FM modulation applied was 
±20MHz around 4. 4022GHz center frequency. The mea- 
sured resonance response appears inverted in the figure 
due to the negative output polarity of the power diode. 
The fact that both bifurcations continue to occur within 
the resonance lineshape (see Fig. I12|) . in spite of the 
short duty cycles that are on the order of ~ ^s, indicates 
that heating processes which have typical time scale on 
the order of s to ms are unlikely to cause these ef- 
fects. Though there are also few who reported a smaller 
time constant of about 3^s @ as a result of local heating 
effects. 

In order to obtain an estimated low limit of local heat- 
ing time scales in our NbN films, we apply a simple hot 
spot heating model to the resonators [23, l28l l29j |. By 
further assuming that the substrate is isothermal and 
that the hot spot is dissipated mainly down into the sub- 
strate rather than along the film [23|, one can evaluate 
the characteristic relaxation time of the hot spot using 
the equation t — Cd/a, where C is the heat capac- 
ity of the superconducting film (per unit volume), d is 
the film thickness, and a is the thermal surface conduc- 
tance between the film and the substrate l^^. Substi- 
tuting for our B2 NbN resonator yields a characteristic 
relaxation time of r ~ 5.4 • 10~®s, where the parameters 
C ~ 2.7 • 10-'''Wcm-3K-i (NbN) d = 3000A (B2 
thickness), and a ~ 1.5Wcm~^K~^ at 4.2K (sapphire 
substrate) have been used. Similar calculation based 
on values given in Ref. yields t ~ 2.1 • O^^s. These 
time scales are of course 2-3 orders of magnitude lower 
than the time scales examined by the FM modulation 
setup. 

Nevertheless intermodulation products observed in 
these resonators l23 | , which generally require fast nonlin- 
ear response for their generation on the order of < 10~"'^''s 
0, further excludes heating effects, though not entirely 
those associated with very fast heating mechanisms sim- 
ilar to hot spots for example. 



E. COLUMNAR STRUCTURE 

It is well known from numerous research works done in 
the past [30I Isij , that NbN films can grow in a granular 
columnar structure under certain deposition conditions. 
Such columnar structure may even promote the growth of 
random built-in Josephson junctions at the grain bound- 
aries of the NbN films. To study the NbN granular struc- 
ture we have sputtered about 2200A thickness NbN film 
on a thin small rectangular sapphire substrate of 0.2 mm 
thickness. The sputtering conditions applied were similar 
to those used in the fabrication of B2 resonator. Follow- 
ing the sputtering process, the thin sapphire was cleaved, 
and a scanning electron microscope (SEM) micrograph 
was taken at the cleavage plane. The SEM micrograph 
in Fig. 1131 which clearly shows the columnar structure of 
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FIG. 11: Frequency sweep time analysis setup. The frequency 
sweep time of the microwave signal generator was FM mod- 
ulated by a saw tooth waveform with frequency /. The re- 
flected power from the resonator was measured by a power 
diode and osciUoscope. 
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FIG. 12: Frequency sweep time measurement. The figure dis- 
plays the resonance measured by Agilent oscilloscope while 
applying a saw tooth FM modulation of frequency 50kHz ( 
Tswecp =20pis) to the signal generator. The left and right bi- 
furcations of the resonance are still apparent in spite of the 
fast rate frequency sweep. Thus indicating that the bifurca- 
tions do not originate from any global heating mechanism. 



the deposited NbN film and its grain boundaries, further 
supports our weak link hypothesis. The typical diameter 
of each NbN column is about 20 nm. 



V. DISCUSSION 

As we argue herein the physical mechanism that may 
potentially account for the nonlinear dynamics of these 
NbN resonators, is the extrinsic weak link mechanism, 
or more specifically Josephson junctions forming at the 
columnar structure boundaries of the NbN films. In ad- 
dition, we show that this hypothesis is qualitatively con- 
sistent with some of our measurement results. 

The measurement results which support the Josephson 
junction hypothesis can be summarized as follows: 



FIG. 13: SEM micrograph displaying a 2200A NbN film de- 
posited on a thin sapphire substrate using similar sputtering 
conditions as B2 resonator. The micrograph exhibits clearly 
the columnar structure of the NbN film and its grain bound- 



(a) The columnar structure of the NbN films and their 
grain boundaries demonstrated in Fig. 1131 which may 
allow the formation of Josephson junctions between the 
NbN grains. 

(b) The onset of bifurcations and nonlinear features 
occurring at considerably low input power levels, about 
2-4 orders of magnitude lower than the onset of Duffing 
nonlinearity in the Nb resonator ( ~ —28 in Fig. |21vcrsus 
^ 10 dBm in Fig. This low power threshold highly 
implies an extrinsic nonlinear origin. Such extrinsic non- 
linearity could be ac driven Josephson junctions, having 
low onset of instability amplitude, or junctions charac- 
terized by a low Josephson critical current, resulting in 
bifurcations in the dynamics of the gauge invariant phase 
difference across the Josephson junctions, at relatively 
low ac current amplitudes [s^- 

(c) The bifurcation dependence on magnetic field as 
was shown in Fig. |U1 where B2 resonance bifurcation van- 
ished as the perpendicular magnetic field was increased 
above some low magnetic field threshold ~ 11.8 mT, 
which is about 3.5 times lower than Hc^ of NbN reported 
for example in On the other hand Josephson junc- 
tions are known for their sensitive dependence on mag- 
netic field. 

(d) The observation of a very small number of bifurca- 
tions in the resonance lineshape, (two in Figs. [3 E and 
three in Fig. |SJ), compared to the large number of built in 
Josephson junctions expected according to the columnar 
structure model, may be explained qualitatively by one 
of the following possible scenarios. The coupling between 
the Josephson junctions and the global resonator acti- 
vates only a small number of Josephson junctions char- 
acterized by some limiting properties such as Josephson 
resonance frequency falling within the resonance band 
of the global resonator, or Josephson onset of instabil- 
ity amplitude being lower than the driving oscillation 



8 



amplitude inside the resonator. Whereas the other sce- 
nario assumes a synchronization mechanism acting on the 
junctions, thus causing them to lock and evolve with the 
same frequency. Such synchronization mechanism may 
be similar to the mechanism recently analyzed in Ref. 
|8,l| , where synchronization between nonlinear oscillators 
is achieved through nonlinear frequency pulling and re- 
active coupling. 

(e) The rapid frequency sweep experiments, which ex- 
cluded most heating effects (except possible very fast lo- 
cal heating mechanisms), further supports the Josephson 
junction hypothesis, since dynamic states of ac biased 
Josephson junctions may involve very little dissipation 
(vanishing dc voltage). 

Moreover it was shown in several experimental works 
published recently 0, |3^ j36„ , examining the reso- 
nance lineshape of a radio frequency resonators coupled 
to a superconducting ring containing a single Josephson 
junction, that interesting nonlinear dynamics can develop 
in the RF circuit resonance response as the external mag- 
netic flux applied through the ring is varied. Among the 
nonlinear dynamics reported therein one can name, bi- 
furcations in the resonance linesha pe, single and opposed 
fold bifurcations (hysteresis loops) |3J|, pinch resonances, 
where the opposed fold resonances appear to touch (pinch 
off) ^3^, and even stochastic jumps ^3^. These experi- 
mental results are shown to be solutions of the nonlinear 
equations of motion for the system. In addition, dif- 
ferent solutions of the nonlinear equation of Josephson 
junction driven by ac field, for various limiting cases, can 
be found in Ref. Whereas the nonlinear dynamics 



of Josephson junctions incorporated in resonant cavities, 
which was analyzed recently in several theoretical and ex- 
perimental studies, mainly for the purpose of developing 
sources of coherent microwave radiation, can be found 
in Refs. j38j, .40^ ,^], and in the cited references 
therein. Although these works do not analyze the effect 
of Josephson junctions on the resonance response of the 
cavity (a single mode is generally assumed), they provide 
a good understanding of the I-V characteristics of these 
arrays. 

In spite of the qualitative arguments presented here, 
a large part of the underlying physics of these nonlinear 
resonators is still not fully understood, and further theo- 
retical work may be needed in order to account for most 
of the experimental data reported here. 
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